Alignment of Rods and Partition of Integers 
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We study dynamical ordering of rods. In this process, rod alignment via pairwise interactions 
competes with diffusive wiggling. Under strong diffusion, the system is disordered, but at weak 
diffusion, the system is ordered. We present an exact steady-state solution for the nonlinear and 
nonlocal kinetic theory of this process. We find the Fourier transform as a function of the order 
parameter, and show that Fourier modes decay exponentially with the wave number. We also 
obtain the order parameter in terms of the diffusion constant. This solution is obtained using 
iterated partitions of the integer numbers. 
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I. INTRODUCTION 

Phase transitions from an isotropic, disordered state 
to a nematic, ordered state are fundamental in equilib- 
rium statistical physics [J] . They occur in liquid crystals 
0, Q , complex fluids , and are closely related to phase 
synchronization 0, 0,13 ■ The standard approach for de- 
scribing order-disorder phase transitions is based on a 
Hamiltonian description. However, it may not necessar- 
ily apply in nonequilibrium settings such as granular sys- 
tems jjj 0, IT(il | , where a non-equilibrium approach, that 
is based on kinetic interaction rules, may be required. 

We study dynamical alignment of rods and show that 
the corresponding kinetic theory framework is appropri- 
ate for describing the underlying order-disorder phase 
transition. This approach conveniently allows character- 
ization of statistical properties such as the distribution 
of rod orientation. 

We consider the mean-field theory for alignment of in- 
teracting polar rods, introduced recently as a model \H\ 
for ordering of microtubules by molecular motors in a 
planar geometry 0, 0, 0, hj,hahj. In this model, 
rods become aligned by pairwise interactions, and they 
may also wiggle in a diffusive fashion, the latter process 
being governed by a white noise. This system reaches a 
steady-state where alignment by interactions is balanced 
by diffusion. With strong diffusion, the system is dis- 
ordered, but with weak diffusion, the system becomes 
ordered. 

Our main result is an exact steady-state solution for 
this kinetic theory. The Fourier modes obey a closed set 
of coupled nonlinear equations. By iterating these equa- 
tions, it is possible to relate all Fourier modes to the 
lowest mode, which in turn is equivalent to the order pa- 
rameter. Similarly, the order parameter itself is shown 
to obey a closed equation and this allows one to express 
the distribution of rod orientations in terms of the diffu- 
sion constant. While in principle, the complete solution 
involves an infinite number of Fourier modes, in prac- 
tice, it is sufficient to consider only a moderate number 



of modes because the Fourier modes decay exponentially 
with the wave number. 

The exact solution is ultimately related to iterated par- 
titions of the integers numbers. This feature is generic: 
it applies to general alignment rules, as well as to arbi- 
trary alignment rates. We also find that depending on 
the alignment rate, the system may or may not undergo 
a phase transition. 

This paper is organized as follows. In Sect. II, we de- 
scribe the rod alignment model. In Sect. Ill, we perform 
the Fourier transform of the governing equation, which is 
then solved in Sect. IV. We discuss nearly perfect align- 
ment in Sect. V, and briefly mention several generaliza- 
tions in Sect. VI. We conclude in Section VII. 



II. THE ROD ALIGNMENT MODEL 



In the rod alignment model, introduced by Aranson 
and Tsimring there is an infinite number of iden- 
tical polar rods, each with an orientation — ir < 9 < tt. 
The rods become aligned via irreversible pairwise inter- 
actions. As a result of the interaction between two rods 
with orientations 9± and 82 , both rods acquire the average 
orientation as follows (see also figure 1): 
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Here, and throughout this study, addition and subtrac- 
tion are implicitly taken modulo 2ir. 




FIG. 1: Illustrating of the alignment process. 

There is also randomness in the form of a white noise: 
each rod wiggles in a diffusive fashion, and this process 
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is characterized by the diffusion constant D. Specifically, 
in addition to the alignment process the orientation 
of a rod changes according to d9/dt = r\ with 7/ an mi- 
correlated white noise (i](t)r}(t')) — 2D6(t — t'). 

Let P(9, t) be the probability distribution function of 
rods with orientation 9 at time t. It is normalized to one, 
fl d6P(9,t) = 1. This distribution function satisfies 
the integro-differential master equation 



dP 
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The first term on the right-hand-side describes diffusion. 
The integral accounts for gain of rods with orientation 
9 as a result of alignment of two rods with an orienta- 
tion difference of <fi, while the negative term accounts for 
loss due to alignment. Without loss of generality, the 
alignment rate is set to 1/2 so that the loss rate equals 
one. This kinetic theory generalizes the Maxwell model 
of inelastic collisions in a granular gas as all rods interact 
with each other and additionally, the alignment rate is 
independent of the relative orientation 
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III. THE FOURIER TRANSFORM 

The governing master equation (0) is nonlinear and 
nonlocal. Its convolution structure suggests using the 
Fourier transform 



P k = (e- lk0 ) = / d9e- lke P(9) 



(3) 



The zeroth mode equals one, Pq = 1, because of the nor- 
malization, and also, Pk = Pl k . The angular distribution 
can be expressed as a Fourier series 
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Since the alignment process Q is invariant with respect 
to an overall rotation 9 —* 9 + ip with an arbitrary phase 
< f < 2tt, if P(9) and Pk are solutions of J3J), then so 
are P(9 + tp) and P k e lkip , respectively. 

The order parameter R = \{e l8 )\ = |-P-i|, with the 
bounds < R < 1 probes the state of the system. A 
vanishing order parameter indicates an isotropic, disor- 
dered state, while a positive order parameter reflects a 
nematic, ordered state. 

We focus on the steady state. Substituting the Fourier 
series (@J into the master equation J2J), setting dP/dt = 0, 
and integrating over <f), we find that the Fourier transform 
satisfies a set of coupled nonlinear equations 



[l + Dk 2 )P k = A-jPiPj 



(5) 



i+j=k 



with Aa = J- r d^e^l" 1 = S^^. The coefficients 



A q are symmetric, A q = A^ q , and explicitly, 
1 



A q = <f0 



kl =2,4,. 

.(-i)^iR M = M, 



(0) 



In J5j, when k = 0, the sum contains only a single term, 
P = P 2 and indeed, P = 1. 

Since Pq = 1 is known, the identical terms A-kPoPk 
and AkPkPo in the steady-state equation (JSJ) are linear in 
Pfc, and thus, we move them to the left-hand side. Then 
Eq. JSJ becomes 

(1 + Dk 2 - 2A k ) P k = J2 A *-i R l P i ■ ( 7 ) 

i+j=k 

Next, we absorb the prefactor on the left-hand side into 
the kernel on the right-hand side by introducing the 
rescaled kernel 
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(8) 



With this transformation, we arrive at a compact, and 
simpler-to- analyze, equation for the steady-state Fourier 
transform 



Pk— £ Gi,j Pi Pj 



(9) 



i+j=k 



The kernel Gij couples the ith and the jth Fourier 
modes. It has the following properties: 



Gij — Gj t i, 
Gij = G-i t -j, 



Gi tj =0, for |i-j|=2,4, 



(10a) 
(10b) 
(10c) 



The governing equation (|5J is nonlinear and moreover, 
for odd fc, the sum contains an infinite number of terms. 
Despite this, it is still possible to solve this equation an- 
alytically! 



IV. EXACT SOLUTION 

First, we notice that Eq. © can be iterated once, lead- 
ing to a sum of products of three Fourier modes 



P k 
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(11) 



Clearly, this procedure can be repeated any number of 
times, leading to a sum over products of any given num- 
ber of Fourier modes. 

Using this iteration procedure for all but the lowest 
Fourier mode, it is possible to express P k as an ex- 
plicit function of P±i only. For instance, when k — 2, 
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Eq. © is simply P 2 — Gi.iPj 2 and when k — 4 one has 
Pi = G 2j2 P|. Thus, the fourth Fourier mode is also 
expressed in terms of the lowest mode, P4 = G 2 , 2 Gf ^P-j 4 
[20j. Odd modes, in contrast, involve an infinite sum. For 

the third mode, P 3 = 2Gi i2 PiP 2 + 2G_i. 4 P_iP 4 H , 

and substituting the second and the fourth moments 
above, 

P3 = 2G1.2G1.1Pf + 2G '_1.4G2.2Gf.1P1P._1 + • ■ ■ . 

Generally, the kth Fourier mode can be written as an 
infinite series involving terms of the form 

pk+npn^ = ph( Pl p_ 1 )n = pk R 2n 

with n a positive integer. Recall that Pi = Re lip . With- 
out loss of generality, we can set the phase to zero, tp = 0. 
Then Pi — R and the Fourier modes can be written ex- 
plicitly in terms of the order parameter R 

00 

P k = R k Y,Pk,nR 2n . (12) 

71=0 

Of course, p ,n = P±i,n = S n<0 . Since p k , n = P-k,n it 
suffices to solve for k > 0. 

The coefficients pk, n represent all iterated partitions of 
the integer k as follows 

fc = 1 + 1 + • • • + 1 + 1-1 (13) 

k+n n 

The iterated partition involves a series of integer parti- 
tions k = i + j, subject to the following restrictions: (i) 
i ^ and j + 1 0, (ii) \k\ ^ 1, (iii) G (J ^ 0. Each such 
partition contributes a factor Gij to p k ^ n (figure 2). For 
example, the partition 2 = (1 + 1) gives p 2 .o = G1.1 and 
the two equivalent partitions 3 = (1 + 2) = (1 + (1 + 1)) 
and 3 = (2 + 1) = ((1 + 1) + 1) give p 3 ,o = 2G 1 , 2 G 1 , 1 . We 
list all non-vanishing terms up to the fifth order in R 



P2.0 = G M , (14a) 

P3,o = 2Gi )2 Gi,i, (14b) 

p 3 ,i = 2G-i A G 2 ,iGl >1 , (14c) 

P4,o = G 2a Gli, (14d) 

p 5 ,o = 2Gx.4G2.2Gf,! + 4G2.3Gi.2G?.!. (14e) 



Here, we used the symmetries (jlOall and (I10b|) to consol- 
idate identical terms. We stress that the partitions used 
here differ from traditional integer partitions in that they 
involve negative integers and in that they are iterated 

m 

Substituting the series expansion l|12|) into the steady- 
state equation JUJl and equating the same powers of R, 
the coefficients pk, n satisfy 

Pk,n = ^ _C PUPj,m- (15) 

l-\-m—n i-\-j—k 
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FIG. 2: The iterated integer partitions. Illustrated is the 
partition corresponding to ^3,1 = 2G_i,4G2,2G?.i. Each par- 
tition k = i + j with i 7^ and j 7^ generates a factor Gij. 
The factor 2 accounts for the two equivalent partitions. 

Since the indexes I and m are positive, this is now a 
recursion equation. Starting with po,o = Pi,o = 1 ; an d 
utilizing the symmetry p k ^ n — P-k.n, equation (| 1 5|) is 
solved recursively. This provides a systematic method 
for obtaining the coefficients pk, n - 

The series (|12|) expresses all the Fourier modes in terms 
of the order parameter R. It remains to obtain the order 
parameter R as a function of the diffusion constant D. By 
Substituting the Fourier solution l|12|) into the governing 
equation © and setting k = 1, we find that the order 
parameter itself can be expressed as an infinite series 

00 

R = Y,r n R 2n+1 - (16) 

n=l 

The coefficients r n are given by the very same recursion 
equation Ijl5|l 

r n = ^2 _1_ Pi,lPi,m- ( 17 ) 

l-\-m—n i+J — 1 

The coefficient r n is the counterpart of the coefficients 
Pk, n and it represents iterated partitions of the number 
1 as in (|13fl . The partitions may not involve 0's. Except 
for the very first partition, the numbers ±1 may not be 
repartitioned (Fig. 2). The first few coefficients are 

n = 2G_ 1 , 2 G 1)1 , (18a) 
r 2 =4G- 2 ,3Gi, 2 G? il , (18b) 
r 3 = 4(G_ 2 ,3G_ 1 ,4 + G_ 3 ,4G 1) 2)G 2)2 G? il . (18c) 

Equation i|16|) is a closed equation for the order param- 
eter. Once it is solved, all Fourier modes follow from 112fl . 
This completes the solution for R and P k , and therefore 
for the angle distribution P(9). 

In practice, one can calculate the order parameter R 
as a root of a polynomial of degree N by truncating i|ltj[l • 
For N — 3, substituting © and (jSJ into (|18a|l yields 
ri = -^(1 + D - 4/tt)- 1 (1 + AD)' 1 , and using lO 
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FIG. 3: The order parameter versus the diffusion coefficient. 
The order parameter was obtained by solving polynomials of 
increasing order. 



FIG. 4: The angular distribution for various values of D. 
The angle distribution was obtained from the first 15 Fourier 
modes. 



gives the cubic equation for the order parameter 
4 1 1 



I? 



3tt D c - D 1 + 4D 



with the critical diffusion constant 



D c = - - 1. 

7T 



(19) 



(20) 



For D > D c , there is only the trivial solution R = 0, 
corresponding to an isotropic state where the rods are 
randomly aligned: Pk = <5fc,o and P(9) = Below 
the critical point, D < D c , there is also the nontrivial 
solution (Fig. 3) 



U= x j^(l+4D)(D C -D). 



(21) 



This corresponds to a nematic phase in which the rods are 
partially aligned. Near the transition point, this align- 
ment is weak, but it becomes stronger and stronger as D 
decreases. The result of Eq. (I21|l is approximate — only 
the cubic term in Ijltjl) has been maintained — yet suffi- 
ciently close to the transition point, the corrections to the 
cubic equation H19|l are negligible and R ~ C(D C — D) 13 
with the prefactor C =y (3/4)(16 - 3tt) and the mean- 
field exponent 0=1/2 [lj,|a,lll| * s as y m Ptotically exact. 
As shown in the Appendix, the uniform state is stable 
for D > D c but unstable for D < D c . 

Well below the critical point, we compute the coeffi- 
cients pk, n and r n to higher order from l|15|) and l|17|l. The 
order parameter is then obtained by numerically solving 
(|TT)|l . truncated at the corresponding order. Since the 
Fourier modes decay exponentially with the wave num- 
ber 



(22) 



and since the order parameter obeys < R < 1, a mod- 
erate number of Fourier modes is sufficient to accurately 



compute P(9). The order parameter rapidly converges 
with N. For instance, N = 11 already provided an accu- 
rate value for R (see Fig. 3). We note that at this order, 
it is still possible to calculate the necessary partitions 
manually. 

Once the order parameter is known, the Fourier modes 
are obtained from Ijl2jl . The steady-state distribution 
becomes 



P(0) 



1 1 

h - R cos ( 

27T 7T 



1 



-G 1A R 2 cos (28) (23) 



-Gi 2 Gi iiT cos (39) + 0(R 4 

7T 



In the vicinity of the transition point, the lowest mode 
dominates. As the diffusion coefficient decreases, the an- 
gle distribution becomes sharply peaked around = 0, 
reflecting that the rods are strongly aligned (Fig. 4) . We 
now analyze this nearly perfect alignment in more detail. 



V. NEARLY PERFECT ALIGNMENT 

For very small diffusion constants, the order param- 
eter approaches unity and therefore, the series solution 
$F2$ is no longer useful. Indeed, very close to D = 0, 
the convergence is slow (Fig. 3). Nevertheless, it is still 
possible to treat the problem by performing the scaling 
transformation 



P(9) 



1 



(24) 



The scaling distribution remains normalized, 
J dxf(x) — 1. At the steady-state, it obeys the nonlinear 
integro-differential equation 



d 2 f(x) 
dx 2 



(25) 



This equation was obtained by substituting l|24|) into the 
master equation setting the time derivative to zero, 
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and replacing the integration limits ±7r/V-D with oo. In 
the limit D — > 0, the periodic nature becomes irrelevant, 
and the problem reduces to the randomly forced inelastic 
Maxwell model, for which several properties including 
the moments and the Fourier transform can be obtained 
analytically [Tlll9l|. 

Consider for example the normalized moments 
M n = (0 n )/n\. They can be expressed in terms of the 
moments of the scaling distribution M n — D n l 2 m n 
with m n = J dxx n f(x)/n\. Substituting this defini- 
tion into the master equation H25|) and performing the 
integration, the moments m n satisfy the closed set of 
equations m n = to„_2 + 2~™ 5Zi_ mi m n-h from which 
the moments are found recursively, too = 1, to,2 = 2, 
and to-4 = =y (the odd moments vanish). The or- 
der parameter can be expressed in terms of the mo- 
ments, R = \ (e i6 ) \ = 1 - M 2 +M 4 H and therefore, 

it is given by a Taylor series in powers of the diffusion 
constant 

1 8 

R= \-2D+—D 2 + --- . (26) 

This expression is useful only for very small diffusion con- 
stants, i.e., D — > 0. 

Since the Fourier transform F(q) = dx e lqx f(x) is 
the generating function of the moments, it also satisfies 
a closed equation, (l + q 2 ) F(q) = F 2 (q/2). Solving this 
equation recursively, the Fourier transform equals the in- 
finite product 18] 



F(q) = 1] 



(27) 



The poles closest to the origin govern the large- x asymp- 
totics. The simple poles located at q = ±i imply an 
exponential decay [lii I22L l23| 

|x| 



/(*) ~ Ce 



(28) 



as | a; | — > oo. The residues to these poles yield the pref- 
actor: C = \ exp n~ 1 {2 2n ~ 1 - l)- 1 ] = 1.47919. 



The Fourier modes satisfy a generalization of (0 
Dk 2 P k = \ (Ai-j+Aj-i-An-A^PiPj. (30) 



i+j=k 



The constant A q is now A q = ± jZ n d</>e iq ^ 2 K((f)). 



A- q , but it is not necessar- 



Again Aq = 1 and A q 
ily true that these constants vanish at even indexes as 
was the case for uniform alignment rates. By excluding 
the zero modes i = and j = from the summation 
in (|30|l . and by following the same steps, we recover the 
governing steady-state Eq. (JSJ . The generalized coupling- 
constants JEl are 



Gi. 



1 A: 



An 



A 2 



A 



2j 



2 1 + £>(*+ j)5 



2A 



(31) 



These coupling constants are manifestly symmetric 
Gi.j — Gj,i. We conclude that the series solutions l|12|) 
and l|16f) hold for arbitrary alignment rates. 

By repeating the steps leading to l|20() , we find that the 
critical diffusivity generally equals D c — 2A\ — 1. The 
condition for having a transition is therefore A\ > 1/2. 
For the constant interaction rate A-i = —, and then 
D c = - — 1. For the hard-sphere rate K((f>) = ~\<f>\, then 
A\ = -tt(2 — 7r), but since A\ < 1/2, the system is al- 
ways disordered. Therefore, depending on the alignment 
rate, there may or may not be an ordered phase. 

We note that in the kinetic theory of gases, analytic 
solutions are feasible only for the Maxwell model, where 
the collision rate is velocity independent p3. EM [2(|, 
while the general Boltzmann equation is not analytically 
tractable. Remarkably, the analytic solution presented 
above does not require a constant rate and for example, 
the hard-sphere like collision rate K(<j>) — C\<j>\ can be 
solved analytically. The discrete nature of the Fourier 
spectrum makes this possible. 



VII. CONCLUSIONS 



VI. GENERAL ALIGNMENT RATES 

The above solution can be generalized in a number of 
ways. For instance, it is straightforward to generalize it 
to imperfect alignment processes, where the orientation 
difference is reduced by a fixed factor (analogous to an 
imperfect inelastic collision), in each alignment event. 

Moreover, it is also possible to analyze situations in 
which the alignment process occurs with arbitrary 
alignment rate K (81 — #2). We assume K{4>) = K(—<fi) 
and without loss of generality, impose the normalization, 
XT d<pK(<j)) — 1. The master equation becomes 

- P(0) [ d4>K{4>)P{6 + ct>). (29) 



In conclusion, we studied kinetic theory for alignment 
of rods. At the steady-state, alignment by pairwise in- 
teractions is balanced by the diffusive motion. At large 
diffusivitics, the system is disordered and at low diffusiv- 
ities it becomes ordered. 

The Fourier modes obey a closed equation that can be 
solved by selective iterations. This allows one to express 
all Fourier modes in terms of the lowest mode, or equiv- 
alently, the order parameter. Similarly, one can obtain a 
closed equation for the order parameter itself, and then 
the solution becomes explicit. Since the Fourier modes 
decay exponentially with the wave number, a moderate 
number of terms is sufficient to accurately obtain the full 
angular distribution. 

This kinetic approach fundamentally differs from the 
traditional Hamiltonian approach for describing phase 
transitions from ordered to disordered states. Yet, the 
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characteristics of the phase transition, including in par- 
ticular, the critical exponent, are identical to the mean- 
field theory. The kinetic theory has the virtue that it 
directly yields distribution functions. 

The exact solution, using integer partitions, is very 
general. It applies to generic "inelastic" alignment rules 
and more importantly, to general alignment rates. The 
periodic symmetry of the problem is ultimately respon- 
sible for this because it involves a discrete Fourier spec- 
trum. It may be useful to generalize this approach to 
alignment of rods in three dimensions. 

The nonlocal and nonlinear governing equation, in its 
simplest form ©, ap pea rs in other physical processes 
including aggregation [22], and it may be interesting to 
utilize the integer partition solution method in other con- 
texts. There are other natural extensions of this work 



including investigation of relaxation toward the steady- 
state, studies of spatial correlations in low-dimensional 
systems, and restricted range interactions where multi- 
ple alignment directions may arise [27L l28l| . 

Experimentally, it may be possible to realize this order- 
disorder transition by vibra ting macroscopic granular 
chains or granular rods 9, 101 where vibration causes 
diffusion and inelastic collisions result in alignment. 
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APPENDIX A: STABILITY OF THE UNIFORM 
STATE 

The uniform state P(0) = j- is always a steady-state 
of the master equation (|2J). To check when this state is 
stable, we write P(0,t) = ^ + p(0,t). To linear order, 
the perturbation satisfies 



dp 
dt 



D 



d 2 p 



p(0-cj ) /2)+p{0 + cp/2) 
2tt 



-p. (Al) 



Let us take a periodic perturbation with wave number k 
and growth rate A, that is, p(9, t) oc e lke + xt . Substituting 
this form into Eq. IjAlfl gives the growth rate 



2A k 



l-Dk 2 . 



(A2) 



The growth rate is positive only for the lowest mode, 
k = 1, and hence, stability is governed by the smallest 
wave number k = 1 for which Ai = D c — D. Indeed, 
the uniform state is unstable below the critical diffusion 
constant l|2(J|) . 



